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We derive a formula for the current through an interacting quantum dot coupled to two su- 
percouducting leads, using the non-equilibrium Green's function formalism. It is shown that the 
formula takes an especially simple form, when the Andreev reflections in the junctions are negligible 
because of large Coulomb repulsion in the dot. Our formalism provides a new framework to inves- 
tigate many-body effects of mesoscopic systems in the presence of the superconducting leads. For 
an example, we describe the resonant tunneling of quasiparticles in the presence of strong Coulomb 
repulsion in the dot. Further, we discuss the boson assisted transport. Multiple sharp peaks are 
found in the current-voltage curve, which originate from the singularity of BCS density of states 
and the boson assisted tunneling. 
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Recent advancement of modern nanofabrication tech- 
nique has stimulated a lot of interest in the study of elec- 
tron transport in quantum dot devicesun. In addition 
to the Coulomb blockade and the resonant tunneling, a 
novel non-equilibrium Kondo effect is expected to occur 
in the low temperature regime for an ultra small quan- 
tum dpLor artificial atom which is coupled to two metallic 
leadscTtEl. More recently, transport through a nanopar- 
ticle is studied experimentally, where the particle is con- 
nected to superconducting (Al) leads rather than normal 
metallic leadsliil. The level spacing of a "superconduct- 
ing" nanoparticle can be either larger or smaller than its 
bulk superconducting energy gap, depending on the size 
of the particle and other parameters such as temperature 
or external magnetic field. When the level spacing of the 
nanoparticle is sufficiently large, the superconducting or- 
der in the dot disapp|ears and the particle can be regarded 
as an ordinary atomlia. 

In this paper, we present a formulation of calculating 
the current through an interacting quantum dot with sin- 
gle level coupled to two superconducting leads (Fig. 1). 
The superconducting leads are characterized by the en- 
ergy gaps Al, Afl, and the chemical potentials Hl, fiR, 
where (ii, — Hr = eV, V being the bias voltage across 
the leads. The subscript L(R) represents the left (right) 
lead. Using the non-equilibrium Green's function tech- 
nique, we derive an expression of the steady state current 
due to quasiparticle processes. It is shown that the cur- 
rent can be expressed in terms of the single-particle spec- 
tral density of states of the dot electron, if the Andreev 
reflection (AR) processes are negligible. In fact, the AR 
contributions to the current can be neglected when the 
Coulomb interaction energy (U) in the dot exceeds both 
the superconducting gaps A l , A# and the coupling be- 
tween the dot and the leads Tl,Tr. It is because, on 
this condition, the 2e-charge fluctuation of the dot is en- 
ergetically unfavourable. Note that this will be a typical 
condition of double barrier tunneling structures with suf- 



ficiently small quantum dot (see e.g. W). In this case the 
formulation becomes quite simple, because there is no 
AC component of the current. 

To describe the system, we start with the following 
model Hamiltonian 

H = H L +HR + H D {4;d a } 

+ E v k (c{ a d a + dUkv) , (i) 

ka€L,R 

where Hl, H-r and Hd represent left, right BCS super- 
conducting leads and the quatum dot, respectively. The 
last term in the Eq. corresponds to the coupling be- 
tween the dot and each of the leads. d\(d a ) creates (de- 
stroys) an electron with spin a in the dot, and ct (cfc CT ) 
creates (destroys) an electron with the momentum k and 
spin er in either the left (L) or the right (R) lead. In 
general, TId includes intra-dot electron-electron interac- 
tions, and coupling to an environment like a bosonic bath 
through the electron-boson interactions!! 

As noted above, when U » Al : r,Tl,r the Andreev 
reflection processes are prohibited, and DC current flows 
by the resonant tunneling of quasiparticles. The current 
flowing out of the lead a (a = L, R) is written by 

where G^ ka (uj) is the Fourier transform of the Keldysh 
Green's function G^ ka (t) = i (c\ a (0)d a (t)^ . By using 

the Dyson's equation in the Keldysh's formalism, the cur- 
rent due to quasiparticles is given by 

^ = i e y " r % mm - 

4- 9 <>)(g;h-g>))], (3) 

where G<(t) = i (dl(0)d a (t)), and G r a and G% are the 
corresponding retarded and advanced Green's functions, 
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respectively. The Green's functions denoted by g are 
the unperturbed Green's functions of the superconduct- 
ing lead. In deriving the Eq.(|^) from the Eq.(||), we have 
neglected terms containing the off-diagonal Green's func- 
tions of dot electron in the Nambu representation because 
those terms should involve Andreev reflection processes, 
which is negligible in our case. It should be noted that 
the "full" Green's functions of the dot represented by G a 
may include the Andreev reflections, but contribution of 
those can be neglected when the charging energy in the 
dot is sufficiently large. Therefore the Eq. (||) can be 
adopted as the formula for the current in the system. 
Using the relations 
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2m [ulf a (E k )8(u - E k ) 
+v 2 k {l- f a (Ek))S{u + E k )] , (4a) 
-2m [u 2 k (l-f a {E k ))5(u;-E k )] 
+v 2 k f a (E k )6(u + E k )] , (4b) 



where u\ 
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2 ^ E k -ti a )> 



represents the (normal) single particle energy of the lead, 
E k = \i a ± ^/ (e k — fi a ) 2 + A|, and f a is the Fermi dis- 
tribution function of the lead a, one find 

x[/ a ( U )((T»-G;( W )) + G<( W )]. (5) 



r a ( a (u! - Ha), with Tq = 

e k ) being the constant characteriz- 



Here Y s a (uj) 

ing coupling between the dot and the lead a which is 
assumed to be independent of the energy, and Ca( w ) be- 
ing the dimensionless BCS density of states of the lead 
a: Cq( w ) — M/v^ 2 — if \uj\ > A Q , and otherwise 
0. It is straightforward to generallize the Eq. (||) to the 
system with a quantum dot containing multi- levels. In 
that case, rf and the Green's functions are replaced by 
M x M matrices, M being the number of the levels, and 
then I a is given by the trace of the M x M matrix in the 
expression ([|). In the steady state, the current is con- 
served, that is Il+ Ir = 0, and thus one can symmetrize 
the formula for the current as 

I = Il = -Ir = \ (h - Ir) 

= j £ /^{(rfM/xM-rfiHMo;)) (6) 
x (g;h - g») + (rf M - rf M) g<m} , 

with the condition of the current conservation (i^ = 
-Ir), 

fdv (rfH + r|H)G<H 
= -/dw (rfM&M + rlH/flM) (7) 

x(G^H-G-H). 



The current formula (^|) is reduced to the simple form 
/= 1?E/ d^r s (u;){f L (uj)-f R (uj)}p a (cj), (8) 

where V s H = ff^ffi^ and p„(u) = -llm G^H, 
by using the approximation 



alto (rf( w )-ii(«)) g<m 



= — I du> 



rfM-r|M 
rfH + r|( w ) 

x (rf (u)f L (u) + T s R (u)f R (u)) (G» - G a a (u)) 



(9) 



This approximation is made by multiplying the factor 



rg(^)-r|(a>) 



on the integrands of both sides of the Eq. (0), 

and seems to be reasonable since 1"? is a smooth 

function of u except for the points uj — [i a ± A a and its 
absolute value is not larger than one. Further, this ap- 
proximation becomes exact in the small coupling limit, 
T <C fcsT, because both G< and G r — G a have sharp 
(5-function-like peaks at resonant level of the dot. 

The Eq. (|J) is the main result of this paper. Note that 
the spectral function p a {^) should be calculated in the 
presence of the leads, and it includes resonant tunneling, 
spin flips, inelastic scattering, etc. For the normal metal- 
lic leads, rf (w) = T Q , and thu^the Eq. (||) is reduced to 
Meir and Wingreen's formulate. That is, the Eq. (||) is 
a generallization of Meir and Wingreen's formula to the 
case of the superconducting leads. 

As an example of the quantum dot, let us consider an 
Anderson impurity. The Hamiltonian of the dot is 



H D = ^s^d^ + Un^n i7 



(10) 



where n a = dl.d a is the number operator of the dot elec- 
tron with spin a. If one neglects higher order effects such 
as the Kondo-like correlations, the Green's function can 
be simply evaluated by the method of the equation 
tion and truncation of high order Green's function! 
To begin with, let us introduce 2x2 Nambu representa- 
tion of the (retarded) Green's function 



G ff (w) =<7 ff 



it >, 



(11) 



where 7 CT 
yields 



d„ 







The equation of motion for G CT 







G a (u) 



(12) 



(2) 



To », 
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where cr z = 



Cka 



-(2) 
la 



_ quantum dot of the Eq. 



n- a d a 



and 1 is the unit matrix. This equation can- 



not be solved exactly due to the last term of the R.H.S. 
In order to get an approximate solution of G CT , we con- 
sider the equation of motion for the Green's function 

< 7. 



(2) 



tions which appear on the equation of motion. Adopting 
this procedure and taking U — > oo limit one can find that 



G„(u) = [G CT M] n 



1 - (n_«r> 



w -e CT + *(rfH + r|H)' 



(13) 



which leads to the resonant-tunneling-like expression of 
the current 

1 = x E (1 - {n - a)) J dLJ [h{uj) ~ Muj)] 
rf(u,)r£H 



( w - £CT ) 2 + (rfH + r|( w ))^ 



(14) 



The factor (1 — (n- a )) arises from strong Coulomb repul- 
sion. Similar expression to the Eq. (|l4] ) has been obtained 
by Yeyati et aZ£3, where the dot Hamiltonian is simply 
simulated by a single non-interacting level, instead of the 
Eq. (|l0|). Resonance broadening of the spectral function 
is included in the Eq. (|l4|). The effect of resonance broad- 
eningiin the I — V characteristic has been discussed in the 
Ref.ta. It is notable that for a small resonance broaden- 
ing r Q <C A a , the current via single level directly reflects 
the BCS singularity in the density of states of the super- 
conducting lead via the function T a (u) (see Fig. 2(a)). 

While the Eq. is valid for A Q > T a , Kondo- 

like correlation becomes important on the other limit, 
A a «C T Q . In order to describe the Kondo effect properly, 
more elaboraterCalculation, such as non-crossing approx- 
imation (NCA)13, is required in calculating the spectral 
function of the dot electron. For normal metallic leads, 
NCA has been used to investigate anomalies of the I —V 
characteristic, arising from the Kondo resonances in finite 
bias voltagestj. The formalism given in this paper can ad- 
dress the problem of the non-equilibrium Kondo effect in 
the presence of the superconducting leads. 

Another interesting example is the boson-assisted 
transport. Consider a quantum dot in which the elec- 
trons interact with bosonic modes w q with the coupling 
strength g q . The Hamiltonian of such a quantum dot is 
given by 



Hd = e a d\d a + Un^n l 

a 

+ ^2?iuj q b\b q + n^ g q (b q +b\), 



(15) 



where n — Y)^ n a and b q (b q ) creates (destroys) a bo- 
son. For the full Hamiltonian Ti of the Eq. (m) with the 



with 4> — i X 
Hamiltonianc 



Lg q /hu) q )(b\ 



dp), 



a canonical transformation 
b q ) gives the new effective 



= H l + Hr+Hd+Ht, (16a) 



■fl and truncate higher order Green's func- where 



Hp = e a n<j + Un^ni + y] Uugb^bg, (16b) 



and 



H T = Yl V *( c L d " ei * + e- i<t '4c ka y (16c) 

ko-£L,R 

The level position and the Coulomb energy in the dot 
are renormalized by the electron-boson interaction as 
ea = e a - £ g .g 2 /^, U = U - 2J2 q g 2 /Jiu; q . The 
phase factors e ±l< ^ which appear in the tunneling term 
are related to the boson assisted tunneling. With the 
transformed Hamiltonian 7Y, the current formula (j^) is 
modified as 



^^{/iM-ZsHJ^H (i7) 



where p a is the spectral function of the composite parti- 
cle denoted by d a = d a e 1 ^ : 



p a (u) = — Im -C d a , dt > 

7T 



In general, exact calculation of p a (ui) is not possible. 
However, in the limit of weak coupling between the dot 
and leads, that is for r = Tl +Tr <C fcsT, the leads can 
be neglected and the spectral function can be evaluated 
exactly for that case. We consider the Einstein model, in 
which all bosons have a same energy hojo and are coupled 
to the electrons in the quantum dot with strength gq_ .In 
this model the spectral function is given by (see e.g.EI) 

oo 

p a (u,)=e-«( 2N °+V J2 / ; (2«[iVo(^o + l)] 1/2 ) 

/ — — OO 

x e m "°/ 2 {(1 - 5(u - e a - Ihcoo) (18) 

+ (n-o-)5(uj - e a -U - lhuj a )} , 

where N a = l/(e /3?iWQ - 1) and a = gl/(fiuj ) 2 are the bo- 
son occupation number and the dimensionless electron- 
boson coupling constant, respectively. The renormalized 
quantum dot level and the Coulomb repulsion can be 
written as e a = e a — cthojQ, U = U — 2ahujQ. h(x) is the 
modified Bessel function. Applying the current formula 
( |l7| ) with ([l8]), one can get I — V characteristics in the 
presence of the electron-boson coupling. 
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Fig. 2 shows the I — V characteristics for different val- 
ues of a. Here we consider the spin degenerate case 
(e-f = ei = e). In the absence of the electron-boson 
coupling (Fig. 2(a)), the I — V curve consists of two sharp 
peaks at eV = 2(\e\ + A) and at eV = 2(e + U + A). 
These peaks reflect the singularity in the BCS density 
of states of the leads. This feature, which makes better 
spectral resolution than the case of the normal metal- 
lic leads, have been verified experimentallyO. Note 
that there is no thermal broadening in contrast to the 
case of the normal metallic leads. For normal metal- 
lic leads, there are two thermally broadend steps at 
eV = 2|e|,2(e + U). For a ^ 0, the main peaks are 
located at eV = 2(\e\ + A),2(e + U + A). Further, we 
obtain additional side peaks at eV = 2(\e + 17ujo\ + A) 
and at eV — 2(e + U + ITiloq + A) with positive integer 
I related to boson emission. Boson absorption processes 
are negligible because ksT <C Tiuiq is taken here. For 
weak electron-boson coupling, only single boson emis- 
sion process contributes to the I — V curve (Fig2.(b)). 
As the coupling becomes stronger, multiple side peaks 
appear (Fig2.(c)). For normal metallic leads, side steps 
arise from the boson assisted tunneling, but those are not 
clearly shown because of thermal broadening. 

In conclusion, we have derived a formula for the current 
through an interacting quantum dot coupled to two suer- 
conducting leads, when the Andreev reflecction processes 
are negligible because of the strong Coulomb repulsion 
in the dot. This formula provides a new framework to 
study interaction effects in the transport through meso- 
scopic systems coupled to superconductors. For example, 
quasiparticle resonant tunneling has been described by an 
equation of motion method in the limit of large Coulomb 
repulsion in the dot. Further, current-voltage charac- 
tersistics of boson assisted tunneling have been studied. 
Splitting of BCS singularities, which arises from the bo- 
son assisted tunneling, has been shown in the I — V char- 
acteristics. 
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FIG. 1. Schematic diagram of a quantum dot connected to two superconducting leads. The superconducting leads are 
characterized by the chemical potentials hl,^r and the energy gaps Al, Ar. An eletron current / flows from left to right if 
fiL — fJ-R = eV > 0. 

FIG. 2. Current-voltage characteristics (a) without, and (b)(c) with boson assisted tunneling. Here Al = A_h = A, 
r_L = r_R = r/2 and fiL = eV/2, fiR = — eV/2 are considered. Other parameters used here are ksT — 0.1, e = e-f = ej_ = —2, 
U — 8, Tiujq = 0.6 in the unit of A. Full lines represent the results for the superconducting leads, while dashed lines being 
corresponding results for the normal metallic leads. 
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